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Abstract 



In this work we study the multilayer adsorption of polyatomic species on 
homogeneous and heterogeneous bivariate surfaces. A new approximate an- 
alytic isotherm is obtained and validated by comparing with Monte Carlo 
simulation. Then, we use the well-known Brunauer-Emmet-Teller's (BET) 
approach to analyze these isotherms and to estimate the monolayer volume, 
v m . The results show that the value of the v m obtained in this way depends 
strongly on adsorbate size and surface topography. In all cases, we find that 
the use of the BET equation leads to an underestimate of the true monolayer 
capacity. 



1 Introduction 



The theoretical description of multilayer adsorption is a long-standing im- 
portant problem in surface science that does not have a general solution. 
[U [2] Mainly this is due to the fact that the structure of the different layers 
differs from that in contact with the solid surface (first layer). At high cov- 
erage (multilayer region), it is expected that the adsorption process is well 
described by the slab theory of Frenkel, Halsey and Hill. 0, HI E] In this 
approach it is assumed that the higher layers retain the structure of the bulk 
liquid, and only its free energy changes gradually as one goes away from the 
solid surface. On the other hand, at low coverage, it is more appropriate to 
use the Brunauer-Emmet-Teller's (BET) isotherm, [6] where the crystal-like 
structure of the surface is considered. In the BET theory it is assumed that 
the molecules are localized in sites and that the adsorption in the first layer 
is different from the remaining ones. 

The BET equation is one of the most widely used isotherms. The ap- 
proach discards such things as the polyatomic character of the adsorbate, the 
interaction between the admolecules and the surface heterogeneities. Then, 
with the purpose of including a more complex situation, numerous gener- 
alizations of the BET theory have been proposed. 0, [7] Nevertheless, the 
simplicity of the BET isotherm has made it very popular for practical pur- 
poses. In fact, by fitting an experimental isotherm with the BET equation, 
in many cases it is possible to predict the monolayer volume (or monolayer 
capacity) of the solid surface with an error not bigger than the 20%. [T] This 
surprising result is attributed to compensation effects arisen as consequence 
of having carried out many approximations, [I] 

By means of numerical experiments, Walker and Zettlemoyer [8] have 
shown that the conventional BET theory predicts a monolayer volume smaller 
than the real value, when heterogeneous adsorption isotherms are analyzed. 
A similar result was obtained by Cortes and Araya, [9] when considering a 
Gaussian distribution of adsorption energies. The authors showed that the 
estimated values of the monolayer volume from the BET equation diminishes 
with increasing degree of surface heterogeneity (the width of the distribution 
of adsorption energies). 

More recently, Nikitas [10] has arrived to similar conclusions by consid- 
ering both, surface heterogeneity and polyatomic character of the adsorbate. 
In ref [TU] , by using an extension of the Flory-Huggins polymer solution the- 
ory, [7j the multilayer adsorption of polyatomic species was studied over a 
random heterogeneous surface. By particular cases analysis, the author con- 
cludes that one can obtain an underestimation of the true monolayer capacity 
of the order of 25%, when the adsorbate occupies more than one lattice site. 
This underestimation is bigger, if an heterogeneous surface is considered. 

In this work, we study how the monolayer volume predicted by BET equa- 
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tion differs from its real value when considering both the adsorbate size and 
the surface topography, i. e. the space distribution of the adsorption energies 
over the solid surface. In particular, we consider the multilayer adsorption of 
polyatomic species on one- dimensional (ID) and two-dimensional (2D) ho- 
mogeneous and heterogeneous bivariate surfaces. In each case, approximate 
analytic isotherms are built and validated by comparing with Monte Carlo 
simulation. Then, we estimate the monolayer volume, by analyzing these 
isotherms with the conventional BET theory. 

The paper is structured as follows. In section 2, we present the lattice- 
gas model. Next, in sections 3 and 4, the multilayer adsorption isotherms 
for homogeneous and heterogeneous surfaces are obtained. The dependence 
of the monolayer volume on the surface topography and the adsorbate size 
is presented in section 5. Finally, conclusions are drawn in section 6. 

2 The Lattice-Gas Model 

A simple lattice-gas model to describe the multilayer adsorption of poly- 
atomic molecules on homogeneous surface has been recently proposed. [TT1 
[T2"] The surface is modeled by a regular lattice of M sites, all with the same 
adsorption energy e, and the adsorbate is represented by /c-mers (linear parti- 
cles that have k identical units). A fc-mer adsorbed on the surface occupies k 
sites of the lattice with an energy ke and can arrange in many configurations. 
This property is called adsorption with multisite occupancy. On the other 
hand, for higher layers, the adsorption of a k-mer is exactly onto an already 
absorbed one, with an adsorption energy of kU. Thus, the monolayer struc- 
ture reproduces in the remaining layers. This phenomenon is called pseu- 
domorphism and, for example, is observed experimentally in the adsorption 
of straight chain saturated hydrocarbon molecules. [13] Finally, following 
the spirit of the BET theory, no lateral interactions are considered and only 
interactions among the layers are introduced. Figure 1 in ref [TT] shows a 
snapshot representing this lattice-gas model. 

We modify the model to consider the adsorption on a heterogeneous sub- 
strate: now the adsorption energy depends on each site i of the surface. 
Then, the Hamiltonian of the system is 

M 

H=(N-N m )kU + J2^i, (1) 

i=i 

where N is the total number of /c-mers, <7j the occupation variable which can 
take the values if the corresponding site is empty or 1 if the site is occupied, 
and N m is the number of fc-mers on the surface (monolayer), 

1 M 

iV m = rE^- (2) 

i=l 



3 



The Hamiltonian in eq ([T]) can be rewritten as 

M 

H = NkU + y £a i (e i -U). (3) 

i=l 

In the following sections, we will look for analytic solutions of the model in 
ID and 2D, for homogeneous and heterogeneous surfaces. 



3 Multilayer Adsorption on Homogeneous Sur- 
faces 

For the previous model, a simple analytic expression of the multilayer ad- 
sorption isotherm can be obtained in few particular cases. If the surface is 
homogeneous (e, = e for all i) and k = 1, it is easy to demonstrate [H [21 [TT] 
that 

a = 1 °P/P0 (A) 

(l-p/ Po )[l + (c-l)p/poY 1 } 

Here, 9 = kN/M is the total coverage, p is the pressure, po is the saturation 
pressure of the bulk liquid and c is a constant defined as 

c = exp [-13k (e - U)] , (5) 

where (3 = 1/ksT is the inverse temperature (being ks the Boltzmann 
constant and T the absolute temperature), eq (JJJ) is the well-known BET 
isotherm [B] and can be applied to systems in any dimension. 

In the case of k — 2, it is only possible to obtain an exact solution in ID 

Hi/2' 



(1 -p/po] 



P/Po 



(Ac - l)p/p 



(6) 



Equation (ED is the exact dimer isotherm for the multilayer adsorption on 
ID homogeneous surface. As it has been previously shown, [TT] the values of 
the monolayer volume predicted by eqs (jlj) and are different: if a dimer 
isotherm is analyzed, the value of the monolayer capacity arising from using 
BET characterization is smaller than the real one. 

In general, the multilayer isotherm corresponding to the model of eq ([3]) 
cannot be expressed by one equation only. To describe the isotherm, it is 
necessary to give two functions. [TUJ |T2] Let us assume that an analytical 
expression of the monolayer adsorption isotherm is known, exp (fi/3) = A, be- 
ing n the chemical potential and A (the fugacity) a function of the monolayer 
coverage 9 m = kN m /M. Then, the following equations can be deduced [T2] : 

P_ = A(fl m ) 
Po c + A (6» m ) 
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Figure 1: Exact ID isotherms for k 
indicated. 



1, 2 and 10, and two values of c as 



and 







0„ 



[i+p/poY 

Equations (JTj) and (jSJ) constitute the multilayer adsorption isotherm. By 
using the surface coverage (0 < 6 m < 1) as a parameter, we can calculate the 
relative pressure from eq ([7]). Then, the values of 9 m and p/po are introduced 
in eq (jHJ) and the total coverage is obtained. 

Following the previous scheme, it is possible to obtain the exact multilayer 
isotherm for /c-mers in ID homogeneous surfaces. [12] We start from the exact 
monolayer isotherm [IT 



A 



Or 



(fc-1) 

k 



Or 



fc-1 



k(i-e m y 



(9) 



Then, eq (jZJ) can be written as 



P_ 

Po 



U 



k 



k-l 



kc(l-0r 



Or 



(k-l) o 
k Ur > 



fc-1 ' 



(10) 



Equations f TTOT) and (IE]) represent the exact solution of the ID model. In 
particular, for k — 1 and k — 2, it is possible to solve these equations 
to obtain single expressions of the multilayer isotherms, eqs (jlj) and (jSJ), 
respectively. Figure [T] shows the exact ID isotherms for k — 1, 2 and 10, and 
two values of c. 

Also, the previous scheme can be used to obtain an accurate approxima- 
tion for multilayer adsorption on 2D substrates accounting multisite occu- 
pancy. In this case, we use the semi-empirical (SE) monolayer adsorption 
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Figure 2: Comparison between analytic and simulated adsorption isotherms 
for 2D lattices and two different values of c as indicated, (a) k = 2 and 6, 
(b) k = 10. In all cases, we have used t = 10 5 MCSs. 



isotherm [15], [16] 



A 



26 „ 



(fc-1) 



(fc-l)0 m r 



2(fc-D fl 



(fc-l)(l-0 m ) 



where 7 is the connectivity of the lattice. It has been shown that eq fill I) is 
a very good approximation for representing multisite-occupancy adsorption 
in the monolayer regime. [T5l [T6] Then, by using eqs (J7J) and ( TTTT) we obtain 



2^,1 



(fc-l)0m 



1 _ 2(fc-l) fl 
1 7 fe 



(fe-i)(i-e m ) 



7 A;c (1 - 9 m ) k + 2^ E 



1 - tt^e 



(fe-i)e n 



■yk U " 



(fc-i)(i-e m ) • 



(12) 

Note that, for 7 = 2, the SE isotherm is identical to the eq Qj. Therefore, 
eqs (JHJ) and (I12p represent the general solution of the problem of multilayer 
adsorption in homogeneous surfaces with multisite occupancy: for 7 = 2 
(ID) this isotherm is exact, but is approximate for 7 > 2. In addition, in 
any dimension, the exact isotherm for k = 1 (BET equation) can also be 
obtained from eqs (jHJ) and (fT2l) . but always with 7 = 2 (these equations with 
k — 1 and 7 > 2 are erroneous, because do not provide the BET isotherm). 

In order to test the 2D approximation, we have compared the analytic 
multilayer isotherm with results of Monte Carlo (MC) simulation. The algo- 
rithm used is described in ref [12] . Here, the equilibrium state is reproduced 
after discarding a number t of MC steps (MCSs). Then, the mean value of 
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the total coverage is obtained as, 



9 



k(N) 



(13) 



M 



where the average (N) is calculated over another t successive MCSs (the total 
number of MCSs is 2t). The computational simulations were developed for a 
square lattice (7 = 4) of linear size L (M = L x L). For each value of k, we 
choose L = 20k. For these lattice sizes (proportional to k), we have verified 
that finite-size effects are negligible. 

Figure EH shows a comparison between the analytic isotherm [given by 
eqs (IE]) and ( fT2l) ] and the MC results, for k = 2 and 6 and two values of c. 
As we can see, the agreement is very good for the parameters used in the 
figure. On the other hand, the accuracy of the analytic isotherm diminishes 
as k increases. Figure [2b shows this effect for k = 10. Also, in this figure, 
we can appreciate that the difference between the analytic and the numerical 
isotherm diminishes as c is increased. 

4 Multilayer Adsorption on Heterogeneous Sur- 
faces 

In the previous section, we have obtained the multilayer isotherm from the 
monolayer isotherm. It is possible demonstrate that in general, the formalism 
allows to establish this connection only if (1) pseudomorphism is present and 
(2) no lateral interactions between the molecules in the multilayer regime are 
considered. In fact, eqs (171) and (El) still hold if the particles in the monolayer 
interact among them and with the solid surface. Although we could use 
this formalism to determine the multilayer adsorption isotherm for a given 
surface heterogeneity (for which it would be necessary to have an appropriate 
monolayer isotherm), we have chosen to use a different strategy. 

We start here from the integral representation of the adsorption multilayer 
isotherm [2] 



where 0\ oc (e) represents the local adsorption multilayer isotherm correspond- 
ing to an adsorptive site of energy e and x( £ ) is the adsorptive energy dis- 
tribution which characterizes the surface heterogeneity (as before, the total 
and the local coverage depend on p and T) . 

It should be noticed that eq ( Tl4l) is strictly and generally valid only for 
noninteracting monomers (k = 1), which is a quite unrealistic case. If ad- 
sorbed particles occupy more than one site (multisite occupancy) or interact 
with each other, then the local coverage at a point with a given adsorptive 




(14) 
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energy depends on the local coverage on neighbor points with different ad- 
sorptive energies and, in general, eq (TH]) should be replaced by a much more 
complex expression. [TT], ITS] 

Nevertheless, in some situations eq (114j) represents a good approximation 
of the adsorption isotherm (see below). For a lattice-gas model of fc-mers, 
we can generalize this equation as 

6 = J20ioc(E s ). (15) 

s 

In the last equation, the sum extends over all possible configurations of a 
single k-mer in an empty lattice, and E s is the adsorption energy of each 
one of them. Note that the values of E s depend, among other things, on the 
energy distribution x( £ )i the surface topography and the number k. 

In following sections, we will study the multilayer adsorption on ID and 
2D heterogeneous surfaces. As local isotherm, we will use eqs (jHJ) and ( fT2l) . 
Then, we will compare the multilayer adsorption isotherm obtained by using 
eq ( TIT)]) and the calculated with MC simulation. 

4.1 Adsorption on ID Heterogeneous Surfaces 

As we said in the introduction, the heterogeneity is modeled by two kinds of 
sites (bivariate surface): strong sites with adsorption energy e\ and weak sites 
with adsorption energy 62 {s\ < 62)- In ID, where the surface is represented 
by a chain of sites with periodic boundary conditions, these sites form patches 
of size I (I = 1,2,3, ■ ■ •), which are spatially distributed (topography) in a 
deterministic alternate way. 

The number of possible configurations of a single fc-mer in an empty 
lattice is M. However, due to periodicity, eq ( JTBTl has only 21 terms (with 
many of them having the same adsorption energy). Then, the multilayer 
isotherm is approximated as 

1 21 

= ^£0ioc(Q). (16) 

Each term corresponds to an effective value of c given by 

Ci = exp [—(3 (Ei — kU)} , (17) 

where Ei is the adsorption energy. This value of c$ can also be expressed as 
function of c\ and C2, the values of c for homogeneous surfaces given by eq <$E§ 
and whose adsorption energies are E\ and 62, respectively. If the i-th term in 
eq ( TLB"]) corresponds to a k-mer with ki units located over strong sites and k 2 
units located over weak sites, then the adsorption energy is Ei = kiei + k 2 £2, 
and 




S 



As mentioned previously, we use eqs © and (THZI) with 7 = 2 as local 
isotherm. Note that eq (fTUj) is a sum of local isotherms with different values of 
c, but at the same relative pressure. Then, in most of the cases it is necessary 
to be careful: although, for each local isotherm the surface coverage should 
be used as a parameter, it is not possible to use this as common parameter. 
In fact, eq ((7|) shows that for a fixed value of p/po, the surface coverage 
depends on c, m (cj). 

Now, we analyze two simple cases. On one hand, eq fTTU]) is exact for 
k = 1 and can be obtained as the semisum of two BET isotherms, 



e 



p/po 



Cl 



C'2 



2 (1 - p/po) { [1 + ( Cl - 1) p/p ] [1 + (c 2 - 1) p/p ] 



(19) 



In this case, eq ( I19p does not depend on / and, consequently, the multilayer 
adsorption isotherm is the same for all topography. 

On the other hand, eq ( TTBT) has three different terms for k = 2, being each 
one of them a dimer isotherm eq ([6]) with a particular value of c. Thus, for 
k = 2, the multilayer adsorption isotherm is 
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+ T 




1 



P/PO 



1 - p/po 



nl/2" 



p/Po) 



1 + (4ci - l)p/p 
1 - p/p 



i + (V 



(1 -p/p ) 



'2 - 1J p/po 

1 - p/po 



1 + (4c 2 - l)p/p 



1/2" 



1/2' 



+ 



+ 



(20) 



The first [third] term in the RHS of eq f[2"Ul) represents the adsorption within 
a strong [weak] patch, on a pair of sites (1,1) [(2,2)], with c\ [c 2 ]. There 
are (/ — 1) configurations of this for each patch. The remaining term of 
eq (l20l) corresponds to a dimer isotherm with adsorption energy E = S\ + e 2 
(c = y/ciOi) . There are only two configurations with this energy for pairs of 
sites (1,2) or (2,1). Contrary to eq (1191) . eq (T2l7|) depends on / and the dimer 
isotherm sees the topography. 

For / = 1, the adsorption energy of a dimer is E = E\ + 62 for all configu- 
ration and, consequently, eq (1201) is exact. In general, for / > 1 and c\ 7^ C2, 
this equation is approximate. Then, to determine the range of validity of 
this equation, we compare the analytic isotherm with MC results. Figure 
shows the dimer isotherm for patches of size I — 2, c\ — 10 2 and different 
values of c 2 . As we can see, for c 2 > 10 the analytic isotherms agree very well 
with the MC data. However, for smaller values of c 2 , the differences between 
theoretical and numerical data begin to be significant. This happens because 
eq (1201) has been built assuming that the three different pairs of sites are filled 
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Figure 3: Comparison between approximate analytic and MC isotherms for 
dimers and different ID heterogeneous surfaces: (a) I = 2, (b) / = 3, (c) 
I = 4 and (d) 1 = 5. In all cases we have used lattices of size L = 1200 and 
t = 10 5 MCSs. 

simultaneously and independently However, for c\ 3> c 2 , the real process 
occurs in 3 stages: (i) the pairs of sites (1,1) are covered; (ii) the pairs (2,2) 
begin to be filled and (in) the multilayer is formed. Note that in the first 
stage all the pair of sites (1,2) and (2,1) are removed. For this regime, a 
better approximation can be obtained by a semisum of two isotherms with 

ci y c 2 . 

When I = 3, Figure [3b, the agreement between the analytic isotherms 
and the MC data is very good for all values of c 2 . In this case, for c\ c 2 
the first stage does not eliminate all the pairs of sites (1,2) and (2,1), because 
each dimer occupies only two sites in the strong patches. For this reason, 
the range of validity of eq ( 1201) is wider than in the previous case. Now, if 
/ = 4 or I = 5, we see in Figures [Hfc and d that the behaviors are similar 
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Figure 4: Function Avsl/k for the adsorption in a ID heterogeneous surface 
with ci = 10 2 , C2 = 1 and different values of k as indicated, (a) Function 
A calculated from eq (|T6|) . The inset shows the dependence of A on c 2 for 
k = 10 and I = 10. b) Function A calculated from eq ([231) . The inset shows 
the dependence of A on C2 for k = 10 and / = 5. 

to those observed for I = 2 or I = 3, respectively. In general, for even Z, 
the first stage eliminates almost completely the pairs of sites (1,2) and (2,1), 
while this does not happen for odd I. Finally, when I — > oo, the fraction of 
pair (1,2) and (2,1) goes to zero and eq (120"]) is exact. This limit corresponds 
to the called large patches topography (LPT), where the surface is assumed 
to be a collection of homogeneous patches, large enough to neglect border 
effects between neighbor patches with different adsorption energies. 

In general, if k ^> I (with k > 1), the multilayer adsorption isotherm can 
be represented by a single homogeneous isotherm 

= e loc (^F 2 ). (21) 

On the other hand, for a LPT where k -C I, the isotherm is 

= ~01oc (Ci) + ^loc (c 2 ) • (22) 

The details of the topography are relevant only when k ~ I. In this case, all 
terms in eq (TL6l) are important. Nevertheless, it is also interesting to obtain 
a simpler expression of the multilayer isotherm given by 

9 = (^r) 6loc (Cl) + G) 6loc (v/ ^ } + \w) 6loc (ca) ■ (23) 

Equation ( l23l) captures the extreme behaviors eqs ( 1221) and ( 12T1) . and it ap- 
proximates the MC isotherm as well as eq (1161) . To verify this statement, we 
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Figure 5: Comparison between approximate analytic and MC isotherms for 
dimers and different 2D heterogeneous surface: (a) I — 2, (b) I = 3, (c) / = 4 
and (d) / = 5. In all cases, we have used lattices of size L = 240 and t = 10 5 
MCSs. 



calculate the integral 



A 



9 



MC 



d(p/po), (24) 



which allows to quantify the difference between the MC isotherm, ^mc, and 
the analytic isotherm, 8, given by either the eq ffT6]) or the new approach 
eq (|2"3"1) . For practical purposes, we have chosen a range of relative pressure 
of — 0.3 to calculate the integral eq f[2"4"j) . MC simulation were carried out 
for lattice sizes of L = 20k with a number of t = 10 5 MCSs. 

Figure shows the function A calculated from eq ( fT6l) for c\ = 10 2 , 
C2 = 1 and different values of k and /. As in the case of dimers, the difference 
between the analytic and the MC isotherms increases when / is approximately 
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Figure 6: Function Avsl/k for the adsorption in a 2D heterogeneous surface 
with ci = 10 2 , C2 = 1 and different values of k as indicated, (a) Function A 
calculated from eq (fl6|) . The inset shows the dependence of A with C2 for 
k = 6 and / = 6. b) Function A calculated with eq (I23p . The inset shows 
the dependence of A with c% for k = 6 and / = 3. 

a multiple of k. However, for large patches, i.e. I > 3k, this difference 
becomes smaller. The inset shows, for a particular case (k = 10 and / = 10), 
how A diminishes as C2 is increased. On the other hand, in Figure Hb we 
can see the function A calculated from eq (I2"3l . For k > 2, the first peak 
is higher than the one shown in Figure Hh and is located in a value of l/k 
between 0.5 and 0.7. Nevertheless, the oscillations attenuate quickly as the 
parameter l/k is increased. As before, but now for k = 10 and I = 5, the 
inset shows how A diminishes as Ci increases. 

Finally, the analysis of Figure H] indicates that, instead of eq ffTUj) which 
has many terms, eq ( l23l) can be used as a more simple expression to approach 
the ID multilayer adsorption isotherm. 



4.2 Adsorption on 2D Heterogeneous Surfaces 

As in the homogeneous case, we represent the 2D surface by a square lattice 
with fully periodic boundary conditions. Strong and weak sites are spatially 
distributed in square patches of size I forming a chessboard. Now, the total 
number of configurations of a single fc-mer is 2M. However, as before, only 
21 terms are necessary to describe the multilayer isotherm. The explanation 
is quite simple: the available energies of a fc-mer that it is forced to move in 
any direction of the lattice (row or column) are the same that in ID. Then, 
eq ( |T6l) continues being valid in 2D, where the local isotherm is given by 
eqs (H2D and ©, with 7 = 4. 
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Figure 7: BET plots for ID analytic isotherms of /c-mers with v m = 1, c = 10 
and different values of k as indicated. Dotted lines correspond to linear fits 
of the data in the range 0.05 to 0.25. 



We begin analyzing the multilayer isotherm for k = 2 (for k — 1, eq fTTUl) 
continues being valid in 2D). Unfortunately, it is not possible to write a 
simple analytic expression [as eq ( 1201) ] in 2D. Nevertheless, the multilayer 
adsorption isotherm for dimers has the same structure that eq (TSUI) . Namely, 
it is composed by three terms with c\, c<i (both multiplied by (/ — 1)/2Z) and 
c = \/c\C2 (multiplied by 1//). Note that in 2D this function is approximate 
for any value of the parameter I. 

Figures b, c and d show the multilayer isotherm for c\ = 10 2 and 
different values of / and ci- The behavior of these curves is very similar to the 
one observed in ID, but the difference between analytic and MC adsorption 
isotherms (for even and odd values of I) is no longer so important. On the 
other hand, Figures [6^ and b show the dependence of the function A on l/k, 
where A was calculated by using eq ffTBl) and eq f[2U]) . respectively. As in the 
case of 2D homogeneous surfaces, the analytic isotherm does not fit very well 
the MC data for k > 6. For this reason, Figures EH and b show the function 
A up to k = 6 only. 

In addition, we have shown that just by using an expression of three terms, 
eq ( l23l) . we can approach very well the multilayer isotherm in ID and 2D for 
the adsorption on heterogeneous surfaces. In the next section, we will use 
this approximation and MC simulations to study how the topography affects 
the determination of monolayer volume predicted by the BET equation. 
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Figure 8: Results of the BET plots for the adsorption in ID homogeneous 
surfaces. Dependence on k of the fractions (a) v^/v m and (b) c*/c, for three 
different values of c as indicated. 

5 Monolayer Volume 

In this section, we carry out numerical experiments to determine, in different 
adsorption situations, how much the value of the monolayer volume predicted 
by the BET equation differs from its real value, v m . With this purpose, 
analytic and MC isotherms were analyzed as experimental data. In this way, 
we have determined how adsorbate size, energetic heterogeneity and surface 
topography affect the standard determination of the monolayer volume. 

In a typical experiment of adsorption, the adsorbed volume of the gas, 
v, is measured at different pressures and at a given fixed temperature. In 
terms of this quantity, the total coverage is 9 = v/v m . Analyzing an isotherm 
with the BET equation, it is possible to estimate the monolayer volume. We 
rewrite the eq (jlj) as 

-P/Po- (25) 



v(l-p/p ) CV m CVv 

This equation is a linear function of p/po- If we denote with a and b, the 
y-intercept and the slope of this straight line, respectively, we obtain 

(26) 



m a + b 
and 

c* = - + 1. (27) 
a 

The asterisk has been added in order to indicate that the quantities given 
by eqs (1261) and (1271) correspond to the prediction of the BET theory. Then, 
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Figure 9: Results of the BET plots for the adsorption in 2D homogeneous 
surfaces. Dependence on k of the fractions (a) ^/w m and (b) c*/c, for three 
different values of c as indicated. 



by means of a plot (the so-called BET plot) of the experimental data of 

v(i-p°/p ) vs p/Poi we can obtain an estimate of the monolayer volume and 
the parameter c. Nevertheless, in the experiments it is commonly found that 
there are deviations from linearity in the BET plot. In many cases, the 
linear range extends from a relative pressure of 0.05 to 0.35, although there 
are cases where the range is shorter, p] 

Although the BET plot is a very simple and popular protocol, the value 
of the monolayer volume obtained in this way can differ from its real value. 
As mentioned in the introduction, in an interesting numerical experiment, 
[H] Walker and Zettlemoyer analyzed a BET plot of an analytic isotherm 
composed by two BET-like contributions (a isotherm similar to eq (1191) for 
a LPT), each one with different values of v m and c. The authors concluded 
that the application of the conventional BET equation to this heterogeneous 
isotherm may lead to an underestimate of the true monolayer volume, with 
a c lying between the values for the two type of sites. Later, Cortes and 
Araya j9] have obtained a similar result by averaging the BET equation with 
a Gaussian distribution of adsorption energy. More recently Nikitas, [TU] has 
arrived to similar conclusions by considering both, surface heterogeneity and 
polyatomic character of the adsorbate. 

In the following, we will show that even for adsorption over homogeneous 
surfaces, the polyatomic character of the adsorbate affects significantly the 
predictions of a BET plot. Next, in Section 5.2, by considering bivariate 
surfaces, we will study the combined effect of energetic heterogeneity and 
multisite occupancy. 
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Figure 10: Results of the BET plots for the adsorption in ID heterogeneous 
surfaces with / = 1 and LPT. Dependence on k of the fractions (a) v^/v m 
and (b) c*/c, for three different values of C\ and c 2 (being Ci/c 2 = 2.5) as 
indicated. 

5.1 Homogeneous Surfaces 

We begin analyzing the BET plots of the multilayer adsorption of /c-mers over 
homogeneous surfaces, given by eqs §8§ and (|T2l) and MC data. Although 
in each particular case it is possible to find an optimum range of relative 
pressures, for practical purposes, we have chosen to set this range from 0.05 
to 0.25. Nevertheless, by choosing other ranges (for example, between 0.05 
and 0.35) we obtain similar results. 

In Figure [7| we show the BET plot for ID analytic isotherms with c = 10 
and k — 1, 2 and 10. Note the deviations from linearity in the isotherms for 
k = 2 and 10, which are concave to the pressure axis. The same behavior 
is observed in experimental isotherms and it is attributed to the existence 
of surface heterogeneities. [8] However, as we see in the example shown in 
Figure El these deviations also appear for the multilayer adsorption with 
multisite occupancy on a homogeneous surface. 

On the other hand, as indicated in Figure [3, the obtained value of v ^ for 
k > 1 is smaller than the real one (we set v m = 1), while the opposite effect 
is observed in the estimate of the parameter c. Figures Eh and b show the 
dependence of these quantities on k for different values of c. In all cases, we 
obtain < v m and c* > c, but the differences between the BET predictions 
and the real values are smaller with increasing c. 

Similar results have been obtained in 2D: the BET plots of both analytic 
and MC isotherms show the same curvature as found in ID. Figures [9h and 
b show the results of these 2D BET plots. As we can see, the differences 
between analytic and MC isotherms are significant for small values of c. 
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Figure 11: Results of the BET plots for the adsorption in ID heterogeneous 
surfaces with / = 1 and LPT. Dependence on k of the fractions (a) v^/v m 
and (b) c*/c, for three different values of c\ and c 2 (being ci/c 2 = 10 3 ) as 
indicated. 

However, always < v m and c* > c for k > 1. As in the ID case, the 
monolayer volume predicted by BET is approximately 10-30 per cent smaller 
than the real value. 



5.2 Heterogeneous Surfaces 

In previous work, [HI E] it has been determined that, as heterogeneous ad- 
sorption isotherms of monomers are analyzed, the monolayer volume obtained 
from a BET plot is smaller than the real value. Since in this case k — 1, the 
surface topography does not affect the obtained results. In this section, we 
study the dependence of the monolayer volume on both, adsorbate size and 
surface topography. In particular, we analyze analytic and MC adsorption 
isotherms of fc-mers over bivariate surfaces with I = 1 and LPT. 

Figures [TUk and b show the results of the ID BET plots for three different 
values of C\ and C2, being C\jci = 2.5. In this case, only analytic isotherms 
were studied because they are exact for LPT and the agreement with MC 
data is seen to be remarkably good for I — 1 [except for odd values of / 
(/ > 1), as was previously mentioned]. In all cases we have used c = y/cic 2 
as the reference parameter. As we can see, the curves show that there is 
not a significant difference between both topographies. Only when the quo- 
tient between c\ and c 2 is increased, the space distribution of the adsorption 
energies over the solid surface begins to be important. This is shown in Fig- 
ures [Hh and b, where ci/c 2 = 10 3 . The results of the BET plots for I = 1 
and LPT are very different. For c\ = 10 and c 2 = 0.01, the deviations due 
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Figure 12: Results of the BET plots for the adsorption in 2D heterogeneous 
surfaces with I = 1 and LPT. Open (full) symbols show results obtained from 
the analytic (MC) isotherms. Dependence on k of the fractions (a) v^/v m 
and (b) c*/c, for three different values of c\ and c 2 (being Ci/c 2 = 2.5) as 
indicated. 

to molecule size k are increased in LPT, i. e. the monolayer volume and the 
parameter c obtained from a BET plot are, respectively, smaller and larger 
than the real values (or the reference value). However, most of the curves 
show a compensation effect which is larger for I = 1, and for c\ = 10 3 and 
c 2 = 1. 

Finally, Figures [T2"h and b, and Figures [THa and b show the results of the 
2D BET plots for ci/c 2 = 2.5 and c 1( /c 2 = 10 3 , respectively. In all cases, we 
have analyzed both analytic and MC isotherms. As we can see, the behavior 
is similar to the ID case. Nevertheless, even taking very different values of 
the parameters c\ and c 2 , it is not possible to obtain a complete compensation 
effect. 



6 Conclusions 

In the present paper, an analytic isotherm for the multilayer adsorption of 
polyatomic molecules on different surfaces has been proposed. The formalism 
reproduces the classical BET theory [6] and the recently reported dimer equa- 
tions; [TT] leads to the exact solution for a ID homogeneous substrate; and, 
as is demonstrated from comparison with MC simulation, provides a good 
approximation for ID heterogeneous surfaces. With respect to 2D substrates 
(homogeneous and heterogeneous surfaces), the approach is not exact. How- 
ever, MC data shows that, for molecules of moderate size (not larger than 
k = 6), the analytic isotherm behaves qualitatively similar to the simulation. 
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Figure 13: Results of the BET plots for the adsorption in 2D heterogeneous 
surfaces with / = 1 and LPT. Open (full) symbols show results obtained from 
the analytic (MC) isotherms. Dependence on k of the fractions (a) v^/v m 
and (b) c*/c, for three different values of c\ and c 2 (being ci/c 2 = 10 3 ) as 
indicated. 

In addition, we carry out numerical experiments to determine, in different 
adsorption situations, how much the value of the monolayer volume predicted 
by the BET equation differs from its real value. For this purpose, analytic 
isotherms and MC data were analyzed as experimental data. For ID and 2D 
homogeneous surfaces, the monolayer volume calculated by the BET plots 
is approximately 10-30 per cent smaller than the real value. On the other 
hand, in all cases, the parameter c* is always larger than c. As the multilayer 
adsorption occurs on a bivariate heterogeneous surface, a compensation effect 
is found, with very different values of c. Nevertheless, in any of the considered 
cases, this compensation is not enough to eliminate the decrease caused by 
the molecular size. 
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